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I.  INTRODUCTION 
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In  spite  of  many  years  of  theoretical  and  experimental  effort,  there  is 
still  no  convincing  picture  of  ring  current  dynamics  during  geomagnetic 
storms.  Even  such  mundane  processes  as  charge  exchange  do  not  seem  to  play 
a simple,  clear  role.  For  example,  the  apparent  agreement  between  observed 
ring-current  decay  and  calculated  charge-exchange  lifetimes  [Uemohn,  1961] 
as  found  by  Swisher  and  Frank  [1968]  and  Smith  et  al . [1975]  has  been 
challenged  by  Tinsley  [1976]  and  by  Lyons  and  Evans  [1976],  partly  on  the 
grounds  that  Liemohn's  lifetimes  are  too  long  (in  view  of  recent  results  on 
neutral  hydrogen  density)  and  partly  on  the  grounds  that  pitch-angle 
distributions  of  low-energy  ions  do  not  behave  in  agreement  with  calculation. 
These  authors  conclude  that,  after  a few  hours  of  recovery  phase,  a large 
fraction  of  the  ring  current  is  He+  (which  for  E < 50  keV  has  a longer 
charge-exchange  lifetime  than  a proton  of  the  same  energy).  The  picture  is 
further  complicated  by  the  discovery  of  energetic  precipitating  0+  ions 
apparently  associated  with  the  storm-time  ring  current  [Shelley  et  al . , 1974] 
which  also  have  a long  charge-exchange  lifetime. 

Charge  exchange  is  an  important  source  of  anisotropy  for  ring-current 
protons.  When  this  anisotropy  becomes  large  enough,  it  leads  to  the  generation 
of  ion  electromagnetic  cyclotron  (EMC)  waves  [Cornwal 1 et  al . , 1970]  which  can 
cause  further  proton  losses  and  possibly  make  SAR  arcs  by  heating  electrons 
[Cornwal 1 et  al . , 1971].  In  this  paper  we  explore  the  interplay  between  the 
charge-exchange  process,  the  rate  at  which  energy  and  anisotropy  is  lost 
during  the  wave-emission  process,  and  the  influence  of  heavy  ions  on  wave 
emission  and  absorption.  Four  major  features  emerge: 

(1)  With  the  aid  of  simplifying  (but  realistic)  assumptions,  a convenient 
analytic  expression  for  the  wave  growth  rate  as  a function  of  time,  of  the 
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parameters  of  the  distribution  function,  and  of  the  cold  plasma  density  is 
found,  for  the  case  of  an  initially  isotropic  distribution.  During  recovery 
phase,  as  the  cold  plasma  density  increases,  the  maximum  growth  rate  occurs 
at  a ratio  of  frequency  u>  to  proton  gyrofrequency  £2  of  w/fi  = 0.4  , 
independent  of  time. 

(2)  The  proton  anisotropy  A grows  (while  the  flux  decays)  to  the  point 
where  finite-amplitude  ion-EMC  waves  are  emitted,  if  the  initial  flux  is 
large  enough.  For  some  time  thereafter,  A is  roughly  constant,  because 
there  is  near-balance  between  charge-exchange  growth  of  A and  decrease  of 

A by  wave  emission.  During  this  phase,  the  anisotropy  is  roughly  governed 
by  the  well-known  expression 

R 8*NA2(1+A)  (J) 

where  E^  is  the  resonant  parallel  proton  energy,  B the  earth's  field,  and 
N the  total  plasma  density.  Relation  (1.1)  is  fairly  well  satisfied  ex- 
perimentally [Mill jams  and  Lyons,  1974a].  The  fact  that  the  anisotropy  ceases 
to  evolve  after  a certain  time  is  an  important  consideration  for  comparing 
observations  with  theory.  Eventually  A begins  to  grow  again,  in  order  to 
keep  up  the  wave  growth  rate  in  the  face  of  a decaying  flux  of  ring-current 
protons. 

(3)  The  rate  at  which  protons  lose  energy  to  the  waves  is  governed  by 
the  charge-exchange  loss  rate,  when  the  only  source  of  anisotropy  is  charge 
exchange.  At  first  sight  this  may  seem  somewhat  surprising,  in  view  of  the 
gross  discrepancy  between  charge-exchange  time  scales  and  wave-growth  time 
scales.  The  simplest  statement  [Cornwa 1 1 , 1975b]  of  this  effect  follows  from 


the  single-particle  emission  picture  put  forward  by  Brice  [1964].  A more 


realistic  treatment  is  based  on  the  quasi-linear  moment  equations  of  Cornwal 1 
[1975a].  We  show  that  sufficient  energy  is  generated  to  power  SAR  arcs. 

(4)  If  heavy  ions  (He+,  He++,  0+)  are  initially  isotropic,  charge  exchange 
will  do  little  to  change  this  during  the  first  few  proton  lifetimes  (except 
for  He++  above  ~ 10  keV).  Consequently,  these  ions  act  to  damp  the  proton- 
EMC  waves,  a heavy-ion  fraction  of  0.2  at  beginning  of  recovery  phase  is 
enough  to  prevent  wave  growth  altogether.  This  damping  is  quite  sensitive 
to  the  very- low-energy  (.<  1 keV)  part  of  the  heavy-ion  distribution  function, 
which  remains  unmeasured.  In  this  paper,  we  do  not  consider  time  scales  so 
long  that  heavy  ions  have  acquired  large  anisotropies  by  charge  exchange. 

A clear  idea  of  the  role  of  charge  exchange  in  ring-current  dynamics  can 
only  be  gotten  at  the  expense  of  ignoring  other  potentially  important  effects. 

For  example,  there  are  most  likely  additional  sources  of  both  proton  flux  and 
anisotropy  (e.g.,  radial  diffusion),  and  there  is  experimental  evidence 
[Williams  and  Lyons,  1974a, b]  that  the  plasmapause  is  a boundary  inside  or  outside 
of  which  certain  effects  do  not  act.  This  evidence  is  that  the  pitch-angle 
distributions  for  E < 50  keV  protons  outside  the  plasmapause  do  not  evolve 
into  the  sort  of  anisotropic  distributions  seen  inside;  rather,  they  become 
isotropic  except  for  the  loss  cone  (flat-topped).  This  could  be  because  of 
an  additional  source  (convection)  or  relatively  weak  pitch-angle  diffusion 
mechanism,  such  as  an  electrostatic  instability,  which  acts  only  outside  the 
plasmapause.  Or  it  could  be,  as  Josel.yn  and  Lyons  [1976]  point  out,  that 
E < 50  keV  protons  pick  up  their  anisotropy  as  they  are  irradiated  by  ion-EMC 
waves  generated  off  the  equator  by  higher-energy  protons,  which  are  continually 
anisotropic.  (We  note  that  the  required  cyclotron-resonant  energy  of  the 
higher-energy  off-equatorial  protons  varies  rather  steeply  with  distance  from 
the  equator,  roughly  like  sin'^9  , where  0 is  the  colatitude.)  The  fact 
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that  we  do  not  consider  such  effects  here  does  not  mean  that  we  consider  them 
unimportant;  we  omit  them  only  for  the  sake  of  clarity. 

Finally  there  are  the  questions  of  what  the  initial  anisotropy  of  the 

low-energy  (E  < 50  keV)  ring  current  is,  at  the  beginning  of  recovery  phase, 

and  whether  these  initial  distributions  evolve  consistent  with  charge  exchange 

3 

losses.  The  general  impression  from  the  S data  for  the  December  17-18, 

1971  storm  [Williams  and  Lyons,  1974a,  b]  is  that  it  is  reasonably  isotropic 
inside  the  plasmasphere  (although  not  really  flat)  and  that  the  anisotropy 
grows  some  during  recovery  phase,  except  at  energies  below  = 10  keV.  (This 
latter  effect  has  motivated  Tinsley  [1976]  and  Lyons  and  Evans  [1976]  to 
suggest  that  this  part  of  the  ring  amount  is  largely  He+  .)  However,  there 
is  by  no  means  quantitative  evidence  that  the  ring  current  10-50  keV  evolves 
consistent  with  an  initially  flat  distribution,  with  charge  exchange  as  the 
only  major  effect,  or  that  the  total  flux  decay  rate  of  the  December  17-18 
storm  is  consistent  with  charge  exchange.  We  persist  with  our  idealized 
scenario  of  initially  flat  pitch-angle  distributions  evolving  according  to 
charge  exchange  because;  (a)  To  some  extent,  the  effect  of  pre-existing 
anisotropies  at  the  beginning  of  recovery  phase  can  be  modeled  by  shifting 
the  time  axis  and  renormalizing  the  growth  rates;  (b)  it  is  important  to  know 
something  about  the  theory  of  the  mutual  influence  of  charge  exchange  and 
wave  emission  on  the  pitch-angle  distributions;  and  (c)  whatever  other  effects 
act  on  the  ring  current,  charge  exchange  is  always  there. 
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II.  GROWTH  RATES  AS  DETERMINED  BY  CHARGE  EXCHANGE 
A.  Contributions  from  Protons 

Because  the  magnetospheric  hydrogen  density  is  greatest  at  low  altitudes, 
the  charge-exchange  loss  rate  A(E,L,ct)  is  a function  of  equatorial  pitch 
angle  a , or  of  y = sina  . It  is  reasonably  well-fit  by  a power  law: 


A = A0(E,L)  y~n  (2.1) 

Liemohn  [1961]  first  calculated  A for  magnetospheric  applications,  and 
found  n = 2 . We  agree  with  Tinsley  [1976]  that  Liemohn's  values  of  A 
are  too  small,  because  he  used  too  small  a value  for  the  neutral  hydrogen 
density.  In  fact,  orbi t-averaged  neutral  hydrogen  densities  calculated  some 
time  ago  [Cornwall  et  al.,  1965]  are  reasonably  close  to  Tinsley's  values, 
and  larger  than  Liemohn's.  The  proton  charge-exchange  rates  based  on 
Cornwall  et  al.  are  shown  in  Fig.  1.  A fit  to  Eq.  (2.1)  gives  a value  of  n 
somewhat  less  than  2 (as  already  noted  by  Cornwall  [1975a]),  but  throughout 
this  paper  we  use  n = 2 because  it  results  in  simple  analytic  expressions 
for  the  ion-EMC  growth  rate.  Note  that  from  Fig.  1 the  time  scale  A"1  for 
protons  has  a minimum  value  of  ~ 4 hours  at  L = 3,  E = 3-15  keV.  The  time 
scale  is  even  less  off  the  equator. 

At  t = 0 (beginning  of  recovery  phase)  we  assume,  in  accordance  with 
the  scenario  of  Section  I,  that  the  protons  of  E < 50  keV  near  the  plasma- 
sphere  (which  has  shrunk  to  L = 3 ) have  an  isotropic  distribution  function. 
An  eyeball  fit  to  smoothed-out  data  of  Williams  and  Lyons  is  roughly  ex- 
ponential with  an  e-folding  energy  Eq  of  ~ 15  keV  at  L = 3 . The 
differential  flux  J evolves  in  time  as 


-7- 


J(E,y,L,t)  = J0(E,L)  exp{-A0(E,L)t  y~£) 


(2.2) 


_E/En 

Jo(E,L)  = Jo(L)  e 0 


(2.3) 


Let  us  calculate  the  growth  rate  y of  proton-EMC  waves,  normalized 


to  the  convective  loss  rate  Vq/2£  » where  v^  is  the  group  velocity  and 
l a length  characteristic  of  the  field  line: 


2y£  _ 21  } 4tt  f e 

W I,  . Mr 


G ; i k I 


Me  G P 


x [ dE  dy  6(y-y  ) [ -f  + i f £ - 1 y 

0 L c w J 


(2.4) 


Here  vp  is  the  phase  velocity,  M the  proton  mass,  and 


2 -.1/2  r E n ,3 

E .2 


(2.5) 


with  Effl  = B /8ttN  as  the  magnetic  energy  per  particle.  The  distribution 


function  f is  related  to  J by  f = M2(2E)-1J  . Using  this,  (2.2)  ana 


(2.3)  in  (2.4)  yields 


2yl  . V2tm  1-x  > f dE  , , f , . I 1-x  ] . , 1 

— l— J T o(L)  L l ^ J o J 

F ^ ft 


x exp{-E/E0  - y'o  AQt) 


(2.6) 


where  tc  = Em(l-x)3x'2,  y2  = 1-(EC/E)  . 


As  t increases,  the  anisotropy  factor  (in  square  brackets)  increases, 


while  the  current  decreases  exponentially.  Therefore  2y£Vg  has  a maximum; 
for  finite-amplitude  ion-EMC  waves  to  occur,  this  maximum  must  be  greater 
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than  unity.  The  integral  in  (2.6)  is  manageable  if  we  approximate  XQ  by 
a constant  independent  of  E (but  depending  on  L ) ; Fig.  1 shows  that  this 
is  reasonable  for  E < 30  keV  . We  then  find,  as  shown  in  the  Appendix, 


M , V11  1,  00  f l-x  1 


exp<-A0t  - EC/E0J T 2X0t  K0U)-Q(a,z) 

(2.7) 


2 = 2|  Vfc],/2  , f Ec  ]'«.  2 

' Eo  J ’ l 1 ' 


(2.8) 


Kq(z)  is  a Hankel  function  of  imaginary  argument.  The  closely-related 
function  Q is  non-standard,  but  easy  to  handle.  It  is  described  in  detail 
in  the  Appendix,  where  it  is  shown  that 


Q(a,z)  + Q(a-1 ,z)  = 2K0(Z) 

(2.9) 

Q(0,z)  = 2Kq(z);  Q(1,z)  = Kq(z) 

(2.10) 

Furthermore,  Q is  a positive  decreasing  function  of  a for  0 < a < °°  , 
all  z , thus  bounded  by  0 and  2KQ(z)  . A simple  approximation  to  Q(a,z) 
which  satisfies  (2.9)  and  (2.10)  is 

2K  (z) 

QU.z)  “ • (2.11) 

Although  strictly  speaking  the  derivation  of  (2.11)  is  only  valid  for  large 
z , (2.11)  can  be  used  with  acceptable  accuracy  (-15%)  down  to  z = 0.1  , and 
is  accurate  to  <3%  for  z - 1 . Substitute  (2.11)  into  (2.7)  to  find  the 
norma1 ized  growth  rate 
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2yi  a 
VG 


2tt2  Jq(L)  MM 
N 


exp{ 


V - W 


(2.12) 


The  first  term  in  square  brackets,  proportional  to  t , represents  the  buildup 
of  anisotropy  from  charge  exchange.  The  term  exp(-XQt)  shows  the  decay  of 
the  total  flux.  The  interaction  of  the  two  terms  lends  to  a growth  rate  y 
which  builds  up  and  then  decays. 

It  seems  natural,  by  comparison  with  the  growth-rate  formula  for  a bi- 
Maxwellian  or  for  a distribution  function  of  the  form  y2\(E)  , to  identify 
the  anisotropy  from  (2.12)  as 


A 


eff 


XQt(l+a)  = 


X t 


o 


(2.13) 


However,  (2.13)  is  both  frequency-dependent  and  dependent  on  Em/EQ  . In 
Section  III,  we  use  another  definition  of  the  anisotropy  A , which  is 
aaapted  to  the  quasi-1 inear  moment  equations  [Cornwall , 1975a],  and  which  is 
independent  of  x and  Em/E0  . One  may  remove  the  frequency  dependence  in 
(2.13)  by  evaluating  the  right-hand  side  at  the  frequency  x at  which  the 
growth  rate  is  maximum;  the  corresponding  Ag^  depends  only  slightly  on 
Em/E0  . (This  is  because  x is  not  far  removed  from  that  value  of  x 
which  makes  E£  = Eq  ; by  the  definition  (2.8)  of  z , it  follows  that  A^ 
depends  mostly  on  XQt  ).  Ae^  does  not  quite  agree  with  A as  defined 
in  Section  III,  but  the  maximum  difference  is  < 0.18  over  the  physically 
interesting  range  of  variation  of  other  parameters.  The  reasonable  agree- 
ment of  these  two  definitions  is  a signal  that  the  quasi-1 inear  moment  equations 
more  or  less  fairly  represent  the  physics  of  a ring  current  decaying  by  charge 
exchange  and  emitting  ion-EMC  waves. 
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Me  now  distinguish  two  different  physical  circumstances  under  which 
(2.12)  might  be  used:  (1)  N is  constant  in  time;  (2)  N increases  with 
time.  Case  (1)  is  appropriate  to  protons  well  inside  the  plasmapause,  while 
case  (2)  is  appropriate  to  protons  just  inside  or  outside  the  plasmasphere , 
as  it  refills  and  expands  during  recovery  phase.  For  simplicity  we  suppose 
that  N grows  linearly  with  t , from  an  initial  value  of  zero: 


N = Nt,  N = const. 


(2.14) 


Observe  that,  from  (2.8)  and  the  expression  for  Ec  , z is  constant  in 
time  (at  fixed  x ) for  case  (2),  and  a decreases  like  t”^  . 

In  Fig.  2,  we  plot  the  expression 
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exp{-Aot  - Ec/Eo) 


(2.15) 


for  the  special  case  Em/Eo  * 0.21  (e.g.  L = 3,  EQ  = 15  keV,  N = 10  ; 
The  normalized  growth  rate  2y£v^  is,  by  (2.12), 


= 2^0  {L)  y 

v-  N o 


(2.16) 


■ «(rf  (tt)  (Jo(L)  - ,0'6) ' 


(2.17) 


-3  2 - 1 

where  N is  in  units  of  cm  , Jq ( L ) in  units  of  (keV  cm  sec  ster)  , 

and  l = LRg  , with  Rg  = one  earth  radius.  According  to  Williams  and  Lyons' 

data,  Jq( L ) is  typically  in  the  range  10^-10^  . Fig.  2 is  useful  for 

calculating  the  growth  rate  in  case  (1). 

In  Fig.  3,  we  plot  the  expression  Y(AQt)"^  , useful  for  calculating 

the  growth  rate  in  case  (2),  for  the  special  case  XQt  Em/EQ  = 0.21  (e.g. 


* 3 3 1 

L = 4,  N = 10J  cm- Jday~  ) . In  this  case,  the  normalized  growth  rate  is 


M 2t r2*f«Aft  ~ 

VG  N 0 oL 


(2.18) 


(0o(L)  ^ 10-*)^ 

o 


(2.19) 


Figs.  2 and  3 show  a rather  similar  change  with  time  of  the  maximum  growth 
rate;  the  main  difference  is  that  for  case  (1)  the  maximum  goes  to  increasing- 
ly large  values  of  x , while  for  case  (2)  maximum  growth  occurs  at  x - 0.35  . 

When  numerical  values  for  JQ (L ) fit  to  the  December  17-18  storm  are 

used  in  (2.17)  and  (2.19),  they  reveal  that  2y£Vg1  exceeds  unity  for  a 

period  ranging  up  to  3A’1  , or  three  charge  exchange  lifetimes.  This  is 

roughly  the  period  of  time  in  which  proton  EMC  waves  can  be  emitted  by  an 
initially  flat  distribution,  although  important  modifications  to  the  growth 
rate  from  heavy-ion  damping  and  from  anisotropy  from  wave  emission  will 
shorten  this  time. 

It  is  most  important  to  note  that  the  formulas  and  figures  of  this 
Section  are  not  applicable  once  Zy-Evg1  exceeds  unity.  As  we  discuss  in 
Section  III,  the  anisotropy  becomes  locked  and  ceases  to  evolve  once  ion-EMC 
waves  are  emitted. 

Even  when  the  initial  distribution  is  not  isotropic,  the  formulas  of 
this  Section  can  be  used  for  case  (1)  to  the  extent  that  the  initial  anisotropy 
can  be  fit  to  the  charge  exchange  law  (2.2)  for  some  value  of  t , say 
t = tQ  . It  is  only  necessary  to  shift  the  time  axis  backward  by  an  amount 
t , and  to  remember  that  JqCL ) no  longer  has  the  significance  of  the 
initial  perpendicular  flux  at  the  equator.  In  case  (2),  it  is  necessary  also 
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to  assume  that  the  initial  plasma  density  is  not  zero,  but  rather  NQto  . 
Normally  this  has  little  effect  on  the  final  growth-rate  formulas,  for 
t >>  to  . 


B.  Contributions  from  Heavy  Ions 

Let  us  calculate  the  contribution  to  2y£v^  coming  from  heavy  ions, 
specifically  He+  , within  the  spirit  of  the  same  approximations  that  went 
into  (2.12).  It  is  no  longer  true  that  Xq(L,E)  is  roughly  constant  over 
the  low-energy  range,  a much  better  approximation  is  that  XQ  is  linear  in 
energy.  For  He+  , a suitable  fit  to  XQ  is: 

I Ua \ X (L ) 

V <UE>  - a;(L)E  i x;CL>  - -z&jgy  . (2.20) 


In  (2.20),  Xq(L)  is  the  proton  charge-exchange  rate,  assumed  independent 
of  energy  for  3 < E < 30  keV  . Above  30  keV,  the  proton  charge-exchange 
rate  falls  sharply,  and  XQ(He)  = XQ(p)  at  E » 50  keV  . Eq.  (2.20)  is 
valid  up  to  about  this  energy. 

For  He+  ions,  the  fundamental  rate  formulas  (2.4)-(2.6)  become 
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(2.21) 


(2.22) 


= 0 - E'/E) 


where,  as  before,  x = w/fi  , M the  proton  mass,  and  n is  the  proton 
gyrofrequency.  We  assume  that  the  helium  ions  are  not  present  in  sufficient 
numbers  to  modify  the  phase  or  group  velocity  appreciably. 

Note  that  for  x = 1/4  , E^  is  much  smaller  than  the  corresponding 
value  Ec  for  protons.  For  x = 0.35  , the  frequency  of  maximum  growth  as 
shown  in  Fig.  3,  E(l/Ec  = 0.02  . Since  the  integral  over  E in  (2.21)  runs 
down  to  E^  , the  He+  distribution  function  (~J/E)  is  probed  at  very  low 
energies,  hundreds  of  eV  or  less.  This  is  also  so  if  the  ions  are  He++  , 
but  the  corresponding  ratio  E^/Ec  for  0+  is  not  as  small  by  nearly  an 
order  of  magnitude.  Even  with  a positive  anisotropy,  heavy  ions  yield  a 
negative  contribution  to  y as  long  as  w > ft(ion)  . A positive  heavy-ion 
contribution  to  y is  possible  for  w < ft(ion)  , if  the  heavy-ion  anisotropy 
is  large  enough.  For  He+  , charge  exchange  probably  does  not  contribute 
dominantly  to  this  anisotropy,  since  a typical  He+  charge-exchange  time 
scale  is  -60  hours  at  the  equator  at  L = 3 . (It  is  about  one  day  for 
0+  at  the  same  place.) 

If  we  were  to  calculate  the  heavy-ion  contribution  to  y using 
-E/E1 

j(i°n)^E,L)  _ j^(l)  e 0 , and  (2.20)  for  XQ(ion)  , a formula  similar  to 
(2.12)  would  emerge.  However,  this  formula  would  be  logarithmically  singular 
at  E^  = 0 , reflecting  the  E-1  singularity  in  the  heavy-ion  distribution 
function  f - J/E  . This  is  no  problem  for  protons,  where  Ec  can  never 
become  too  small,  but  the  singularity  in  f must  be  removed  for  the  heavy 
ions.  Let  us  (in  the  absence  of  very- low-energy  measurements)  assume  the 
He+  current  J’  to  be,  at  t = 0 : 


J'(L,E,t=0)  = j;(L)  |r  e 

o 


(2-24) 
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This  Maxwellian  phase-space  distribution  function  may  be  used  either  for 
the  thermal  component  or  for  the  ring-current  component  (in  which  case  pre- 
sumably E^  = Eq  ).  In  analogy  with  the  calculation  of  (2.12),  we  find, 
using  (2.20)-(2.24) 
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(2.25) 


where  the  K^(z)  are  Hankel  functions  of  imaginary  argument.  In  (2.25), 
the  following  parameters  appear; 


Ec 
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z' 


1 1/2 
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1 + Xt 
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(2.26) 


This  expression  is  well-behaved  at  E^  = 0 ; it  is  worth  noting  that  the  heavy- 
ion  contribution  to  the  growth  rate  does  not  decrease  exponentially  at 


First,  consider  the  contribution  of  ring-current  (>  keV)  He+  ions  . 
For  all  practical  purposes,  we  may  set  E^  = 0 in  (2.25),  and  ignore  the 
terms  in  square  brackets;  this  approximation  breaks  down  only  after  many 
proton  charge-exchange  lifetimes  or  for  very  small  x . Then  we  find 
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Then  let  us  take  = Eq  = 15  keV  . The  term  in  square  brackets  is  nearly 
equal  to  unity  for  many  proton  lifetimes.  Finally,  let  us  take  the  He+ 
flux  at  E = Eq  and  at  t = 0 to  be  constant  fraction  of  the  proton  flux 
for  all  L , that  is,  J^L)  = C JQ(L)  • A very  simple  modification  to  the 


1 ! 1-X  1 

growth-rate  formula  (2.15)  emerges:  simply  subtract  ^ C ' ^rom  the 

right-hand  side  of  this  expression  for  y . This  may  be  done  directly  on 
Figs.  2,  3.  Since  5 might  range  from  0.05  to  0.2,  this  can  be  a very 
significant  negative  contribution  to  the  growth  rate.  In  fact,  if  t,  > 0.2 
2y £vl^  will  never  exceed  unity,  and  proton  EMC  waves  will  not  be  generated 


For  0+  , the  charge-exchange  rate  is  also  rather  small,  about  an  order 
of  magnitude  less  than  that  for  protons.  We  may  roughly  incorporate  the 
effect  of  0+  ions  by  taking  5 to  be  the  ratio  (He+  + 0+)/p  . 

For  He++  , the  situation  is  somewhat  different.  Between  10  and  30  keV, 

+*f  + 

the  charge-exchange  rate  for  He  -*■  He  is  comparable  to  the  proton  charge- 
exchange  rate.  The  He++  contribution  is,  crudely  speaking,  given  by  (2.12) 
except  that  inside  the  square  brackets  1-x  is  replaced  by  0.5-x  , and  Ec 
is  found  from  (2.23)  using  i2(He++)/£l  = 0.5  ; the  values  of  a,z  from  (2.8) 
are  correspondingly  modified.  In  contrast  to  He+  , the  He++/p  ratio 
(at  least  from  10  to  30  keV)  is  roughly  constant  in  time.  With  c(He++)  _ 0.05 
to  0.2  , these  ions  contribute  modestly  to  growth  for  x = 0.35  , when 
AQt  > 2 , but  their  main  effect  is  to  build  up  the  He+  population.  This 
can  be  accounted  for  by  letting  c(He+)  be  time-dependent. 

As  for  the  thermal  heavy-ion  contribution,  it  is  negligible,  since 


Ec  >>  Eo  ^or  Practica^  purposes. 
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III.  COUPLED  CHARGE-EXCHANGE  AND  ION-EMC  EMISSION  EFFECTS 


To  understand  this  Section,  the  reader  will  have  to  shift  mental  gears 
from  Section  II,  where  growth  rates  was  calculated  without  considering  the 
effect  of  wave  emission  on  the  distribution  function  f . The  point  is 
that  f will  become  more  anisotropic,  thus  increasing  y according  to 
(2.12)  only  up  to  the  time  where  R = 2y£v^  reaches  unity.  Let  us  call 
this  time  t-|  . After  t1  , finite-amplitude  ion-EMC  waves  are  emitted,  and 
it  is  well-known  that  a necessary  condition  for  wave-particle  equilibrium  is 
R = 1 for  those  frequencies  at  which  finite-amplitude  waves  are  present 
[Cornwal 1 , 1966,  1975a;  Kennel  and  Petschek,  1966].  (Note:  R = 1 is  only 
appropriate  if  there  is  no  heavy-ion  damping  of  the  waves  as  they  pass 
through  the  ring  current.  If  the  heavy-ion  damping  rate  is  _YHe>  R = 1 
should  be  replaced  by  the  larger  value  R = 1 + 2yHe£Vg~^  . This  qualifica- 
tion should  be  kept  in  mind  whenever  we  discuss  the  condition  R = 1 .)  For 
conditions  appropriate  to  the  December  17-18  storm,  Figs.  2 and  3 and  Eqs. 
(2.17)  and  (2.19)  suggest  that  XQt^  can  be  as  small  as  0.2,  which  corresponds 
to  fairly  small  values  of  the  anisotropy  (to  be  defined  below). 

The  constancy  of  R in  time  for  t > t ^ is  achieved  through  a balance 
of  the  rates  of  change  of  the  shape  of  the  distribution  function  in  y , of 
the  cold  plasma  density,  the  decaying  omnidirectional  flux,  etc.  It  is  an 
enormously  complicated  job  to  evaluate  the  precise  y-dependence  of  f which 
(a)  keeps  R = 1 over  a finite  frequency  band  width,  (b)  is  consistent  with 
particle  losses  from  charge-exchange  and  pitch-angle  scattering.  (See 
Etcheto  et  al . [1973]  for  an  attempt  to  solve  a similarly  complicated 
problem.)  We  try  a much  more  modest  program  [Cornwal 1 , 1975a],  in  which  the 
shape  of  the 


anisotropy  A . Now  it  is  no  longer  possible  to  make  R = 1 over  a finite 
frequency  bandwidth;  instead  we  impose  only  the  condition  R = 1 , where  R 
is  the  value  of  R at  the  frequency  x which  maximizes  R . (The  idea  is 
that  the  fastest-growing  waves  largely  determine  the  equilibrium  between 
wave  growth  and  loss.) 

It  is  especially  important  to  realize  that  A is  a dynamical  parameter, 
no  longer  given  by  a simple  formula  like  (2.13).  That  is,  A responds  not 
only  to  charge-exchange  growth  (as  in  (2.13))  but  also  to  loss  of  anisotropy 
by  wave  emission  [Brice,  1964]  and  to  increase  of  anisotropy  from  precipita- 
tion into  the  loss  cone.  Indeed,  A must  respond  very  rapidly  by  wave 
emission  in  order  to  keep  R at  its  equilibrium  value  of  unity;  no  other 
parameter  (e.g.,  EQ,  Em>  XQt)  can  change  at  a rate  comparable  to  the  wave- 
growth  time  scale. 

We  have  noted  before  [Cornwal 1 , 1975a]  that,  of  all  the  quantities  on  which 
R depends,  it  depends  most  sensitively  on  A (at  least  if  A is  less  than  one). 
Thus  the  condition  R = 1 serves  to  determine  A , primarily,  and  the 
constancy  of  R suggests  a similar  near-constancy  for  A . As  a first 
approximation,  anisotropy  growth  from  charge-exchange  balances  anisotropy 
loss  from  wave  emission,  and  the  shape  of  f evolves  much  more  slowly  during 
the  period  of  wave  emission  than  before.  This  approximation  breaks  down 
also  at  later  times,  when  the  flux  has  decayed  so  much  that  greater  and 
greater  anisotropy  is  needed  to  keep  growth  going. 

We  turn  to  the  calculation  of  R . The  first  problem  is  how  to  para- 
metrize the  shape  of  f by  a single  parameter  A , in  order  to  calculate 

R . The  actual  shape  is  unknown,  and  such  a parametrization  will  be  useful 

% 

only  if  similar  results  are  achieved  for  a variety  of  (smooth)  shape  functions. 

2A 

For  a shape  given  by  a power  law  f ~ y or  a bi-Maxwellian  (A  = (T  /T(|  )-l)  , 
numerous  results  are  available  in  the  literature.  Numerical  results  for  R 


are  given  in  Cornwal 1 [1975a]  for  the  bi-Maxwell ian  case.  In  the  present  paper 
we  consider  instead  distribution  functions  whose  pitch-angle  shape  is  given 
by  (2.2),  the  charge-exchange  form,  except  that  X t is  replaced  by  a 
dynamical  parameter  t equivalent  to  the  anisotropy: 

J(E,y,L)  = J0(L,t1)  eT  exp{-E/EQ  - xy"2}  . (3.1) 

Of  course,  at  the  time  t^  when  finite-amplitude  wave  emission  begins, 
x = XQt^  , but  thereafter  x (or  A ) is  determined  by  the  condition  R = 1 . 

How  is  x related  to  the  anisotropy?  For  purposes  of  the  quasi-linear 
moment  equations,  the  anisotropy  must  be  defined  in  terms  of  the  expectation 
values  <NE>  and  <NEX>  of  the  total  energy  and  of  the  perpendicular 
energy.  The  expectation  values  would  be  correctly  defined  by 

<Q>  » | dy  dJ  fQ  . j dEdy  yT(y)  QE1/2  f (3.2) 

for  any  function  Q . The  integrals  are  taken  over  the  first  two  invariants 
y,  J at  fixed  L , and  T(y)  is  the  normalized  bounce  time  as  defined  by 
Schulz  and  Lanzerotti  [1974].  Unfortunately  the  curvature  of  the  field  lines 
obscures  to  some  extent  some  simple  physical  points  we  wish  to  make  later  on, 
so  we  define  expectation  values  as  simply  integrals  over  the  equatorial 
velocity  components  (as  if  the  Earth’s  field  lines  were  straight).  Then  the 
anisotropy  A is  defined  as: 


<NEi>  = <NE> 


A + 1 


(3.3) 


The  relation  between  x and  A is  shown  in  Fig.  4;  note  that  A grows  at 
roughly  twice  the  rate  x does.  A comparison  of  A as  defined  by  (3.2) 
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with  (see  (2.13))  shows  discrepancies  at  the  lowest  values  of  A , but 

these  are  unimportant. 

Given  the  distribution  function  (3.1),  the  calculation  of  R (or 
equivalently  the  maximum  growth  rate)  is  the  same  as  the  calculations  of 
Section  II,  except  that  At  is  replaced  by  t , and  except  that  f no 
longer  decays  in  time  in  quite  the  same  way  as  it  would  if  only  charge 
exchange  were  acting. 

To  explain  this  latter  point,  we  observe  that  if  A is  roughly  constant 
in  time,  protons  must  be  scattered  in  pitch  angle  through  a substantial 
fraction  of  90°  in  one  charge-exchange  lifetime  so  that  anisotropy  loss  from 
pitch-angle  scattering  can  balance  anisotropy  growth  from  charge-exchange. 

In  turn,  this  means  that  protons  do  not  experience  charge-exchange  losses 
appropriate  to  a fixed  pitch  angle,  but  rather  that  all  protons  are  lost  by 
charge-exchange  processes  at  a pitch-angle  averaged  rate  A , defined  by 


A - A 


(3.4) 


A is,  in  fact,  the  loss  rate  which  occurs  in  the  quasi-1 inear  moment- 
equation  for  <N>  or  <NE>  as  we  discuss  below.  Fig.  5 shows  a plot  of 
AA~^  vs.  A . Note  that  A is  roughly  twice  XQ  at  small  A , but  at 
large  A there  are  very  few  protons  at  small  y , and  here  A = Aq  . In 
addition,  protons  are  lost  by  precipitation  into  the  loss  cone  due  to  pitch- 
angle  scattering,  and  energy  is  lost  by  wave  emission.  These  effects  may  be 
incorporated  by  allowing  JQ(l,t^)  and  Eq  into  (3.1)  to  be  functions  of 
time,  whose  time  dependence  is  governed  by  the  quasi-1 inear  moment  equations 
of  the  next  Section  as  well  as  by  the  condition  R = 1 . 
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( 


Figure  5.  Pitch-angle  averaged  charge-exchange  rate  A vs.  anisotropy  A . 
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The  upshot  of  all  this  discussion  is  the  following  formula  for  R : 


R = 2-rr2  J^L.t,)  ^ MAXx{ 


1 -x 
x 


exp[-E  /E  -A(t-t,)] 
1+a  -*  c 0 1 


(3.5) 

where  JQ(L,t1)  is  the  flux  factor  corresponding  to  the  time  t^  at  which 
wave  emission  starts  (i.e.,  R first  attains  the  value  of  unity).  In  (3.5), 
2 and  a are  calculated  as  in  (2.8),  with  x substituted  for  At  (ex- 
cept in  Ec  ).  As  a function  of  A , R shows  the  same  features  as  it  would 
for  a bi -Maxwellian:  a steep  dependence  on  A for  A < 1 . This  is 

illustrated  in  Fig.  6,  where  the  A-dependence  of  R is  isolated  in  the  ex- 
pression 


fMAX 


_ A(t-t, ) 
R e 1 


2tt2  J0(L,ti) 


(3.6) 


for  Em/EQ  ~ 0.21  • The  interested  reader  should  compare  Fig.  6 to  Fig.  1 
of  Cornwall  [1975a]  which  gives  R for  a bi-Maxwel lian;  the  A-dependence  of 
R^  is  quite  similar  for  the  two  cases. 

In  the  next  Section,  we  use  R = 1 as  a dynamic  equation  relating  the 
parameters  which  appear  in  (3.5).  This  equation  is  adjoined  to  the  quasi- 
linear  equations  for  the  distribution  function,  yielding  a set  of  equations 
from  which  the  emitted  wave  energy  can  be  calculated,  as  well  as  the  evolu- 
tion of  A,  <N>  , and  <NE>  . 
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MAX 


IV.  DYNAMICS  ACCORDING  TO  THE  MOMENT  TRANSPORT  EQUATIONS 


In  addition  to  the  equation  for  R = 1 discussed  in  the  last  Section, 
there  is  the  quasi-1 inear  diffusion  equation  expressing  the  rate  of  change 
of  f as  determined  by  charge  exchange,  loss-cone  precipitation,  and  wave- 
induced  diffusion.  We  replace  this  equation  by  a set  of  moment  equations 
[Cornwall , 1975a]  for  the  moments  <N>,  <NE>  , and  <NEX>  . Actually  we  need 
only  consider  equations  for  <NE>  and  <NEX>  , since  <N>  does  not  appear 
in  these  two  moment  equations. 

The  main  objectives  of  this  Section  are  to  determine  the  rate  at  which 
particle  energy  is  lost  to  waves,  and  to  discuss  the  evolution  of  the  aniso- 
tropy during  the  period  of  wave  emission.  It  will  be  very  helpful  to 
consider  a simplified  version  of  the  energy-loss  problem,  following  Brice 
[1964].  He  points  out  that  the  conservation  of  energy  and  momentum  during 
the  elementary  process  of  emission  of  a wave  quantum  by  a single  particle 
lends  to  a relation  between  the  rate  at  which  the  particle  loses  energy 
to  waves  and  the  rate  EiW  at  which  it  loses  perpendicular  energy: 

&W/£iW  = u/n  (4J) 


Assume  that  E/E^  = <NE>/<NEi>  ; then  with  the  aid  of  the  definition  (3.3) 
of  the  anisotropy  A , we  find 


A = 


— E - 
E-E, 


(4.2) 


It  is  simple  to  combine  (4.1)  and  the  time  derivative  of  (4.2)  to  find 


x A, 


W 


A + f 


x 


-1 


(4.3) 
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where  x = w/ft  . Here  Aw  is  the  rate  at  which  anisotropy  is  lost  by 
wave  emission.  Note  that,  since  the  emission  frequency  is  bounded  by 
x < A(l+A)_1  , the  square  brackets  in  (4.3)  are  always  positive.  If  (as 
argued  in  the  last  Section)  A is  roughly  constant,  then  -Aw  must  roughly 
be  equal  to  A^  , the  rate  at  which  A increases  by  charge  exchange.  For 
a distribution  function  of  the  type  (3.1),  we  have 


A = 4*  i = x — =2X 
mCEX  dx  o di  ‘Ao 


(4.4) 


where  we  consulted  Fig.  4 for  dA/dt  . Thus  (4.3)  yields  £W/E  in  terms  of 

XQ,  A , and  x . For  typical  values  of  A and  x , Ew/E  is  15-20%  of 

\Q  ; in  words,  the  particles  lose  energy  to  waves  at  ~15%  of  the  rate  at  which 

they  lose  energy  by  charge  exchange.  We  shall  see  below  that  this  rate  is 

somewhat  larger  when  (4.3)  is  modified  to  account  for  precipitation  losses,  but 

even  this  rate  is  enough  to  power  SAR-arcs.  For  example,  during  recovery 

phase  of  the  December  17-18  storm,  the  total  energy  stored  in  the  ring 

current  is  of  the  order  5 x 10  ergs  cm  (energy  per  unit  area  at  the  iono- 

-2  -1 

sphere),  and  it  is  dissipated  at  a rate  of  = 0.3  ergs  cm  sec  by  charge 

-2  -1 

exchange.  Twenty  percent  of  this  is  0.06  ergs  cm  sec  , while  it  is  estimated 

-2  -1 

by  Williams  et  al . [1976]  that  roughly  0.05  ergs  cm  sec  is  needed  to 
power  the  observed  SAR-arcs  during  this  storm.  The  SAR  arc  begins  to  fade 
rapidly  about  10  hours  (or  2X'1 ) into  recovery  phase. 

Eq.  (4.3)  does  not  take  into  account  a number  of  effects;  to  remedy 
this,  we  turn  to  the  moment-transport  equations.  We  refer  the  reader  to 
Cornwal 1 [1975a]  for  all  details.  The  equation  for  <NEX>  reads 

<NEA>  = -<XNEx>  - 2y  W x_1  . (4.5) 
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W is  essentially  the  energy  density  of  the  waves.  Note  that  there  is  no 
precipitation  loss  term,  since  the  equatorial  value  of  E is  very  small  in 
the  loss  cone. 

The  first  step  is  to  express  all  the  terms  in  (4.5)  in  terms  of  A and 
<NE>  . Since  A = A^"2,  E = Ey2,  <ANEX>  = AQ<NE>  . The  process  is  com- 
pleted by  using  the  definition  (3.3)  of  A , and  we  express  the  result  as  an 
equation  for  2y  W : 


2y  W 


1 

f A + 1 

{ -VNE> ' 

<NE> 


A + f 


(4.6) 


The  idea  is  to  combine  (4.6)  with  the  moment  equation  for  <NE>  , eliminate 
9/3 t <NE>  , and  thus  arrive  at  an  equation  for  2y  W which  contains  no  time 
derivatives.  The  <NE>  moment  equation  is 


<NE>  = -A<NE>  - 2y  W - <A  NE>  (4-7) 

where  we  have  used  (3.4)  to  write  <ANE>  = A<NE>  . The  last  term  in  (4.7) 
represents  precipitation  into  the  loss  cone  produced  by  wave  diffusion.  By 
eliminating  2y  W from  (4.7)  we  find 


<NE>  - <Ar 


(4.8) 


| T <NE>  = d'U  -(A-xA  )<NE>  + 


A 


2 


The  inequality  follows  from  x < A/(A+1)  . Alternatively  we  eliminate 
d/3t  <NE>  from  (4.6)  to  find 


27  W = d-y[  A (A  * Pt  - A ]<NE>  + ^4  <A 

tL  A + 4 0 -i  L A + 4 J F 


NE>  - 


A<NE> 
2 A + | 


(4.10) 

The  analog  of  -£WE~^  in  the  simplified  Brice  model  (4.3)  is  2y  W<NE>  1 . 
It  is  easily  checked  that,  if  the  charge-exchange  and  precipitation  loss 
terms  in  (4.10)  are  dropped,  (4.10)  reduces  to  (4.3)  (with,  of  course, 

A = Aw)  . Furthermore,  the  term  in  square  brackets  in  (4.10)  represents  the 
growth  of  anisotropy  from  charge  exchange.  To  see  this,  substitute 
3/at  <NEi>  = -Aq<NE>  and  3/3t  <NE>  = -A<NE>  in  the  time  derivative  of  the 
definition  (3.3)  of  A , and  come  to 


2 A + | 


A + 1 


A + ! 


(4.11) 


which,  in  turn,  is  easily  shown  to  be  identical  to  (4.4)  with  the  help  of 
(3.1),  (3.3),  and  (3.4).  Since  A = kQEX  + Ay  , (4.10)  can  be  written 


2y  W = d 


■M- 


VNE>  A + 1 

2fA  + 4l2  lA  + ! 


<A  NE> 

P 


(4.12) 


which  is  again  the  simple  Brice  relation  (4.3)  except  for  the  <A  NE> 


term. 


It  only  remains  to  discuss  <Ap  NE>  , the  precipitation-loss  term 
coming  from  pitch-angle  scattering.  We  have  already  argued  that  a particle 
changes  its  pitch  angle  at  about  the  charge-exchange  rate,  so  we  expect 
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<A  NE>  ~ X <NE>  or  possibly  A<NE>  . This  means  that  the  waves  are  in  the 
p o 

weak-diffusion  regime,  because  the  characteristic  pitch-angle  diffusion  time 

-1  — 2 
XQ  is  considerably  smaller  than  the  average  minimum  lifetime  = l^/2aQ 

where  is  the  bounce  time  and  aQ  the  loss-cone  pitch  angle.  However, 
it  is  possible  that  there  is  a transient  regime  in  which  strong  diffusion 
occurs,  if  <NE>  at  the  time  t-j  (when  the  waves  switch  on)  is  quite  large 
compared  to  a critical  value  (discussed  in  Cornwal 1 [1975a])  which  is  roughly 
comparable  to  the  stably-trapped  limit  of  Kennel  and  Petschek  [1966].  This 
transient  regime  ends  when  strong  diffusion  has  reduced  <NE>  to  a value 
comparable  to  the  critical  value  <NE>C  . 

In  order  to  discuss  both  the  strong-diffusion  transient  and  the  weak- 
diffusion  regime  together,  it  is  necessary  to  take  into  account  the  non- 
linear dependence  of  <Xp  NE>  on  2y  W . In  an  earlier  work  [Cornwal 1 , 1975a] 
we  introduced  the  form 


<A„  NE>  = 2y  W<NE>  <NE> 

P |_  C 


_ A.  + I ) x .j.  2v  W 1 

777|  MIN  J • 


(4.13) 


The  term  linear  in  W (found  by  setting  TMJN  = 0 ) yields  the  weak-diffusion 
limit,  while  in  the  strong-diffusion  limit  <Xp  NE>  -*■  <NE>  T^^  . The  critical 
value  <NE>£.  is  approximately 


„ B2  f *p  I dl  A + f 


c 47r  fc  > (A  + 1 ) 


(4.14) 


Use  (4.13)  in  (4.10)  to  find 


2y  W = ^4  IT  Q + (Q2  + 4Tmjn  Y<NE>  )1/21  (4.15) 

2dT„...  A + l > L MIN  c J 
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where 


Q * 

<N£>  - 

<NE>c  + 

'kin 

r A + \ 1 

Y = 

t>- 

L 

' A + 1 J 

A 

0 J 

<NE>  - 


A<NE> 


2(A+1 ) A + | 


(4.16) 


(4.17) 


From  (4.15)  various  limits  can  be  found.  When  TMIN  Y » ) <NE>  - <NE>c|  , 
we  find  that  2y  W is  given  by  the  right-hand  side  of  (4.10),  with  <Xp  NE> 
set  equal  to  zero.  When  <NE>  - <NE>C  » Y , we  find 


2y  W = 


dT 


MIN 


A + 1 


<NE>  - <NE> 


A + ~ 


(4.18) 


This  is  the  limit  which  was  studied  earlier  [Cornwall , 1975a].  It  amounts  to 
saying  that  precipitation  losses  are  much  more  important  than  charge-exchange 
losses,  which  is  true  for  several  special  cases.  The  first  case  is  that  of 
electron  EMC  waves  generated  by  energetic  electrons,  for  which  there  is  no 
charge  exchange.  The  second  case  is  the  transient  regime  mentioned  above  for 
the  proton  ring  current,  when  the  waves  switch  on  in  a very  strong  ring 
current.  The  waves  grow  rapidly,  and  precipitation  losses  reduce  <NE>  - <NE>C 
to  order  of  T^^  Y , at  which  time  all  quantities  begin  to  change  more  or 
less  on  the  charge-exchange  time  scale. 

We  may  estimate  the  size  of  2y  W in  the  weak-diffusion  regime  by 
supposing  <NE>  - <NE>c  - T^^  Y . There  is,  under  these  circumstances,  a 
term  in  (4.15)  of  order  A0(Aq  » P^us  Terms  of  0(Aq<NE>)  . 

Since  a typical  value  of  (A  T|^)”‘  is  about  5,  this  square-root  term  is 
not  an  order  of  magnitude  greater  than  the  terms  of  0(Aq)  , and  the  whole  of 
(4.15)  is,  practically  speaking,  0(Aq<NE>)  . In  other  words,  the  precipitation 


-34- 


loss  terra  <Ap  NE>  in  (4.12)  is  of  order  Aq<NE>  , as  would  be  expected  if 

the  pitch-angle  diffusion  coefficient  is  0(Aq)  • 

As  time  goes  on,  another  sort  of  weak-diffusion  regime  sets  in.  In  this 

regime,  it  is  no  longer  true  that  A is  small  and  that  Y * 0(Aq<NE>)  . 

Rather,  Y = 0(TMIN  ^10<N^>)  » which  means  that  the  anisotropy  is  growing  at 

very  nearly  the  charge  exchange  rate  (cf.  (4.17)  and  (4.11)).  Also, 

|<NE>  - <NE>C[  = °^0TMIN<N^>c^  * and  W = o(Aq<NE>)  . Ultimately,  as  the 

flux  decays,  this  ordering  of  <NE>  - <NE>c  will  be  violated  and  one  finds 
2 — 

2y  W = Aq  TM[N<NE>  . At  this  point  wave  emission  has  effectively  ceased,  if 
it  has  not  already  been  forced  to  cease  because  it  is  no  longer  possible  to 
satisfy  R = 1 . 

The  quantitative  statement  of  these  effects  requires  simultaneous  solution 
of  the  equation  R = 1 (see  Eq.  (3.5))  and  the  9<NE>/3t  Eq.  (4.7)  which  can 
be  written,  using  (4.13)  and  (4.15),  as 


Q + 


(Q2  ♦ 4T„in  Y<NE>c)'/2 


(4.19) 


(4.20) 


Although  some  preliminary  efforts  have  been  made  in  this  direction,  they  do 
not  add  appreciably  to  our  qualitative  understanding  as  described  earlier. 

This  understanding  may  be  summarized  as  follows: 

There  are  three  terms  in  2y  W as  given  in  (4.10).  The  first  (in 
square  brackets)  corresponds  to  the  growth  of  anisotropy  from  charge  exchange. 
The  second  (precipitation-loss)  term  is  estimated  to  be  = Aq<NE>  . The 
third  term  involving  -A  is  smaller  than  the  other  two  terms  until  wave 
growth  is  about  to  end,  when  it  nearly  cancels  the  first  term  (i.e.,  Y in 
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(4.17)  becomes  small).  During  the  main  wave-growth  period  when  A = 0 , 
the  precipitation-loss  term  may  add  significantly  to  2y  W , perhaps  even 
doubling  the  result  coming  from  charge  exchange  alone.  That  is,  the  naive 
Brice  relation  (4.3)  may  only  yield  one-half  the  energy  lost  to  waves,  which 
could  thus  be  20-40%  of  the  energy  lost  by  charge  exchange.  In  view  of  the 
imperfectly  efficient  processes  by  which  wave  energy  is  used  to  heat  electrons, 
this  should  be  approximately  right  to  drive  SAR-arcs. 
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APPENDIX 


The  mathematical  task  is  to  evaluate  Eq.  (8),  for  constant  XQ  . The 
relevant  integral  is 


fCO 

F 


X t E 


£-11 

e-eJ 


A.t  E„ 
0 c 


(A.  1 ) 


The  change  of  variables 


q - ± (E-Ee) 
0 


(A. 2) 


yields 


j n 


(A. 3) 


With  the  aid  of  a standard  formula  for  Hankel  functions  of  imaginary  argu- 


kvU>2) 


bv  dt  t+v_1  exp- 
■ 0 


1 t + £ 
2 1 t 


(A. 4) 
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(A. 5) 


(A. 6) 
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Evidently  Q(0,z)  = 2KQ(z)  , from  (A. 4),  and  Q is  a positive  decreasing 
function  of  a for  positive  a . In  (A. 6),  the  change  of  variables  q -*•  q 
easily  yields,  along  with  (A. 4), 

Qta.z)  + Q(a_1,  z)  = 2Kq(z)  (A. 7) 

from  which  Q(1 »z)  = Kq(z)  . One  may  give  a formal  expansion  of  (A. 6)  as  a 
power  series  in  a (or  in  a~^  ),  with  coefficients  K^Cz)  > but  this  is 
not  very  useful.  By  writing  the  denominator  of  (A. 6)  in  the  form 

dAe'Uq+a)  (A. 8) 

0 


1 

q+a 


some  useful  forms  can  be  derived.  We  quote  only  one  of  them: 


QCa.z)  = 2Kq(z) 


(A. 9) 


Finally,  a very  simple  and  useful  approximation  follows  from  the 
steepest-descent  evaluation  of  (A. 6)  for  large  z . It  is: 


Q(a,z) 


2Kq(z) 

1+a 


(A. 10) 


This  satisfies  (A. 7)  for  all  z . The  approximation  (A. 10)  is  an  overestimate 
of  Q for  0 < a < 1 . It  is  useful  even  for  small  z ; the  maximum  error 
is  about  15%  for  z = 0.1  . 
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LABORATORY  OPERATIONS 


The  Laboratory  Operation*  of  The  Aeroapace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  arid  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 


Aerophysic 8 Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 

fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine 


Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting. superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forma  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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